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ABSTRACT 

In many theoretical scenarios it is expected that intermediate-mass black holes 
(IMBHs, with masses M ~ 10 2-4 M ) reside at the centers of some globular 
clusters. However, observational evidence for their existence is limited. Several 
previous numerical investigations have focused on the impact of an IMBH on the 
cluster dynamics or brightness profile. Here we instead present results from a 
large set of direct N-body simulations including single and binary stars. These 
show that there is a potentially more detectable IMBH signature, namely on 
the variation of the average stellar mass between the center and the half-light 
radius. We find that the existence of an IMBH quenches mass segregation and 
causes the average mass to exhibit only modest radial variation in collisionally 
relaxed star clusters. This differs from when there is no IMBH. To measure 
this observationally requires high resolution imaging at the level of that already 
available from the Hubble Space Telescope (HST) for the cores of a large sample 
of galactic globular clusters. With a modest additional investment of HST time to 
acquire fields around the half-light radius, it will be possible to identify the best 
candidate clusters to harbor an IMBH. This test can be applied only to globulars 
with a half-light relaxation time < 1 Gyr, which is required to guarantee efficient 
energy equipartition due to two-body relaxation. 

Subject headings: stellar dynamics — globular clusters: general - - methods: 
n-body simulations 
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Introduction 



Theoretical work has suggested that some globular cluste rs may harbor intermediate - 



mass black holes (IMBHs; M ~ 10 2-4 M & ) in their centers (e.g.. |Portegies-Zwart et al.ll20041 ). 
If this is indeed the case, there are significant consequences for ultra-luminous X-ray sources, 
gravitational wave em ission from dense st a r clusters, and the dyn amics of globular clusters 
(GCs) in general (see Ivan der Marell 120041 ; iMiller fc Colbertl 12004 fo r an overview). Defini- 
tive evidenc e for IMBHs has, however, been elusive. For example, iGebhardt. Rich, fc Ho 
(120021 . 120051 ) argued for an IMBH in Gl based on the an alysis of HST line - of-sigh t veloc- 
ity data and Keck spectra, but an alternative analysis by iBaumgardt et al.l (J2003aj) points 
out that accep t able dynam ic models without a large central object also fit the observations. 
Gerssen et al.l ( 120021 . 120031 ) argued that the kinematics of M15 seem to slightly favor the 



presence of an IMBH, b ut for this cluster alternative interpretations exist (IBaumgardt et al. 



2003bl ; iDull et al.ll2003l ). More recently, the observed line-of- sight kinematics of Omega Cen 
have also been used to argue for the presence of an IMBH (INoyola. Gebhardt. Bergman 
2008h . 



A more secure identification of an IMBH in a GC can, in principle, be provided by 
also measuring the proper motion of central stars in order to reconstruct their orbits and 
thus firmly establish if a central massive point object is present. Several HST-GO programs 
based on this idea have been approved in past cycles (e.g. GO10474, PI Drukier; GO10401 
& GO10841 PI Chandar; GTO/ACS10335 PI Ford), but to date they have not yielded any 
indisputable detections. The limitation for such studies is the need to carry out multi-year 
observations, thus progress is slow. To maximize the chances of success it is thus of primary 
importance to focus the observations on the candidates most likely to harbor an IMBH. 

Candidate selection is possible if one focuses on the in direct influence of the I MBH 
on the dynamics of its host. Direct N-body simulations by IBaumgardt et al.l (120041 ) and 
Trent i et al.l (l2007bl ) found that the presence of an IMBH acts as a central energy source 
that is able to prevent gravothermal collapse and thus maintain a sizable core to half-mass 
radius ratio throughout the entire life of the GC. The existence of such a large (> 0.1) core 
to half-mass radius ratio in a collision ally relaxed cluster might be due to the presence of 
an IMBH (see also iHeggie et al.l 120071 ). However, the picture becomes more complicated 
when this signature is transferred from the ideal world of N-body simulations, where a 
complete knowledge of the system is available, to real observations, where essentially only 
main seque nce and red g iant branch stars define the light profile of the system. In fact, an 
analysis by iHurleyl (120071 ) cautioned that the difference between mass and light distributions 
can lead to a large observed core to half-light radius ratio for GCs with single stars and 
binaries only. 
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Here, we continue the search for indirect IMBH fingerprints by focusing on the conse- 
quences of the presence of an IMBH on mass segregation. Through direct N-body simulations 
we show that the presence of a large (of order 1% of the total mass), central mass signif- 
icantly inhibits the process of mass segregation, even among only visible main sequence 
stars and giants. To th e best of our knowledge this effect was first briefly mentioned in 



Baumgardt et al.l (120041 ). but left without further quantitative analysis. Quenching of mass 



segregation is present in all of our simulations with an IMBH, independent of the initial con- 
ditions of the cluster, including variations in initial mass function, density profile, strength 
of the galactic tidal interaction, number of particles and initial binary fraction. We find that 
a differential measurement of the average mass between the center and the half-light radius 
is effective in separating star clusters with and without an IMBH, provided that the stellar 
system is at least 5 initial half-mass relaxat ion times old. This me asure is observationally 
feasible with current data (for example, see iDe Marchi et al.l 120071 and references therein), 
and can lead to the selection of a promising set of IMBH host candidates. A direct obser- 
vational application of this approach is left to a companion paper. Here, we focus instead 
on building the theoretical framework for such analysis. In § 2 we describe our numerical 
simulations, in § 3 we discuss our results, and in § 4 we present our conclusions. 



2. Numerical Simulations 



The numerical simulations presented in this paper have been carried out wit h a state-of- 
the-art direct N-body code for st ar cluster dynamics, NBODY6 (jAarsethl 120031 ) . NBODY6 
has been modified as discussed in iTrenti et al.l (l2007bl ) to improve accuracy in the presence 
of an IMBH, and uses regularization of close gravitational encounters without any softening. 
This makes it optimal to follow interactions within the sphere of influence of the IMBH. 



2.1. Units and timescales 



NBODY6 uses the standard system of units of iHeggie fc Mathieul (119861 ) in which G = 
M = —4Et = 1, where G is the gravitational constant, M is the total mass and Et is the 
total energy of the system. In this system of units, the half-mass relaxation time, which 
i s the relevant timescale for mass segregation and energy equipartition is defined as follows 
JSpitzer!ll987i : 



t 



0.138iVr 



rii 



3/2 
h 



In (O.lliV) 
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where N is the number of st ars in the system and is the half-mass radius. In physical 
units t r h can be expressed as ( IDjorgovski Ill993l ): 



L rh 



8.9- 10 V (1M Q 

in(o.niv) x V(^) 



x 



M 



0.5 



lpc 



1.5 



(2) 



where (m*) is the average mass of a star. 



2.2. Initial conditions 



Galactic GCs are made of some of the oldest stars in our galaxy (e.g., see iKrauss fc Chaboyer 



20031 ) and are collisionally relaxed s ystems. Their two-bo dy half-mass relaxation times t r h 



are shorter than their age (e.g., see iHeggie fc Hutl 120031 ). thus their initial conditions are 



largely unknown. However, the evolution on a relaxation timescale is only weakly depen- 
dent on the initial configuration as the system evolves toward a self-similar configuration 
where the density and light profiles are determined primarily by the efficiency of kinetic 



energy production in the core due to gravitational encounters ( jVesperini fc Chernoffl [1994 



Trenti et alil2007aB ). In this paper, we explore a number of different initial configurations, 
varying the initial mass function, the initial density profile, the strength of the galactic tidal 
field, and the fraction of primordial binaries to verify that the evolution of the system is 
indeed inde pendent of t he initial configuration. The initial density profile is always that of a 
single-mass iKingl ( 119661 ) model, but we use a full mass spectrum in the N-body calculations. 
The number of particles is varied from iV = 8192 to N = 32768 to quantify the evolution of 
mass segregation with and without an IMBH. The details of our runs are reported in Table [TJ 



Initial stellar masses were drawn from either a ISalpeterl (119551 ) or iMiller fc Scald (119791 ) 
initial mass function (IMF, hereafter), that is: 



£(m) oc m c 



(3) 



with a = -2.35 and m e [0.2 : 1OO]M for the Salpeter IMF, while for the Miller & Scalo 
IMF the power law slope is the following: a = —1.25 for m e [0.2 : 1]M , a = —2.0 for 
me [I: 2}M Q , a = -2.3 for m G [2 : 1O]M , a = -3.3 for m 6 [10 : 100]M Q . In addition, 
we have also carried out control runs that extend the IMF down to O.1M , as such stars 
exist in GCs, but are not bright enough to be detected in most of the cluster. 

We handled stellar evolution by assuming a turnoff mass M t o . = O.8M , and instanta- 
neously reducing all stars to their final state at the beginning of the simulation. Stars with 
masses O.8M < m < 8.0M^ were assum ed to become white dwarfs and reduced to a final 
mass as prescribed in lHurley et al.l (120001 ). Stars in the ranges 8.O-25.OM and 25. 0-100. OM 
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became neutron stars (unity retention fraction assumed) and black holes, and were reduced 
linearly to 1.3-2.OM and 5-10M Q , respectively. Our model makes the approximation that 
most of the relevant stellar evolution occurs on a timescale shorter than a relaxation time. 
This choice is appropriate to model the dynamics of an old GC on a relaxation timescale 
with only a limited number of particles and is more realistic th an using an unevo lved mass 
spectrum appropriate for young star clusters when N < 30000 ( Trent i et al.ll2008l ). 



In addition, about half o f our runs included primordial binaries, an important compo- 
nent of many GCs (e.g., see iPulone et al.l 120031 ) that can influence the evolution of mass 
segregation. In fact, binaries are on average twice as massive as singles and thus tend to 
segregate in the core of the system (e.g. 



sec 



Heggie et al.ll2006l ). We define the fraction of 



binaries to be 



f b = n b / (n s + n b ) 



(4) 



where n s and n b are the initial number of single stars and binaries, respectively. Thus, a 
run with N = n s + n b = 8192 and f b = 0.1 actually has 8192 + 819 = 9011 objects. As the 
dynamical influenc e of binaries tends to saturate for f b « 0.1 (jVesperini &: Chernofilll994j ; 



Heggie et al.ll2006l ). all our runs with primordial binari es have f b = 0-1, a number similar 



to the observe d binary fract i on of many old GCs (e.g., lAlbrow et al.ll200ll ). Binaries were 



initialized as in lHeggie et al.l (120061 ). i.e. from a flat distribution in binding energy from e r 



to 133e min , with e r 
the cluster. 



(m*)er c (0). Here, cr c (0) is the initial central velocity dispersion of 



To half of the simulations, we added an IMBH with mass Mjmbh ~ 0.01 (~ 1% of the 
entire cluster), the same ratio as a ~ 10 3 M Q black hole would have to a GC of mass 10 5 M Q . 
In some of the simulations (see Tab. [TJ, we increased the mass of the IMBH to study the 
dependence of mass segregation on this parameter. 

All objects in our runs were treated as point masses, thus neglecting stellar evolution 
and collisions as well as any growth of the IMBH due to accretion of tidally disrupted stars. 
These effects have only a minor influence on the late-time dynamics of the cluster, since 
actual collisions are rare after massive stars have evolved, and accretion onto the IMBH is 



minimal (e.g. see iBaumgardt et al.ll2004l ). 



The evolution of the clusters includes the tidal force from the parent galaxy, assuming 
circular orbits with radii such that the tidal cut-off radius is self-consistent with the value of 
the King parameter W used. The galactic tidal field is treated as that due to a point mass, 
and the tidal force acting on each particle is computed using a linear approximation of the 
field. Particles that bec ome unbound are re moved from the system. For full details of the 
tidal field treatment see lTrenti et al.l (j2007al ). 
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For validation p urposes, we also analyzed a few snapsho ts from three run s with N = 
131072 carried out by iBaumgardt and Makino I (120031 ) and by iBaumgardt et al.l (12004 ) with 
and without a central IMBH, kindly made availab le by the authors. These runs include 
full stellar evolution using the iHurley et al.l (120001 ) tracks, but no primordial binaries. In 
addition, accretion of tidally disrupted stars onto the IMBH is included. The initial star 
positions and velocities from these runs were also drawn from a K ing model with Wq = 7.0, 
but instead the mass spectrum was drawn from a iKroupal (120011 ) IMF, with a = —1.3 for 
m G [0.1 : O.5]M and a = —2.3 for m > O.5M . The upper cut-off mass was 15M for 
run 128kk.l, 3OM for run 128kkbh.la, and 1OOM for run 128kkbh.lb. These snapshots 
provide us with a control group against which we can test the validity of our own models, 
and also allow us to probe the extrapolation of our results to higher N. 



3. Results 



3.1. Overall evolution of the star clusters 



A star cluster with single stars only and no central IMBH evolves toward core collapse 
within a few relaxation times. The collapse is eventually halted when the central density is 
high enough to dynamically form binaries. At this stage gravothermal oscillations set in and 
the density profile of the cluster remains self-similar until the final stages of tidal dissolution 
of the system. The existence of either primordial binaries or an IMBH serves as an energ y 
source to counteract the collapse, resulting in a more significant core (jTrenti et al. 2007b ). 
This is confirmed in all models (for example, see Fig. [ J for the evolution of the core and half 
mass radius in our 32k simulations). Here, we use the ICasertano fc Hutl (119851 ) definition of 
the core radius, namely 



z^=i r iPi m i 

L<=i pi m i 



(5) 



where mi is the mass of the ith star, is its from the cluster center of mass, and the density 
Pi around each particle is calculated using the distance to the fifth closest neighbor. 



As expected from previous investigations based on equal mass particles (jTrenti et al. 



2007a]), the density profile of our clusters progresses to a self-similar configuration, which is 
independent of the initial configuration of the stars and the IMF. This independence justifies 
our treatment of stellar evolution at the beginning of the simulations, and provides further 
evidence of the erasure of initial conditions after a few relaxation times. The overall evolution 
of star clusters with and without an IMBH and /or primordial binaries has been discussed 



Trenti et al. 2007a b 



in the literature (IBaumgardt and Makino 



Fregeau &: Rasio 



2003 



2007 



Bau mgardt et all 12004 ; iHeggie et all 12006 



Hurleyll2007l ). Here we focus instead on a novel 
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aspect that has a promising observational signature, namely the evolution of mass segregation 
for runs with an IMBH. 



As our overarching goal is to propose a viable observational test to identify a star clus- 
ter that is likely to harbor an IMBH, we took steps throughout our analysis to replicate 
observational data as closely as possible. We limited our analysis to data projected onto 
two dimensions, and excluded stellar remnants from our calculations of observationally ac- 
cessible quantities. Most of our runs did not have stars with masses below O.2M , because 
these are generally too dim to detect with a high completeness. However, we did perform 
two 16k control runs with masses down to 0.1M Q to ensure that the presence of smaller, 
undetectable stars did not affect the mass segregation. For a proper comparison to the other 
runs, we did exclude the stars with masses between 0.1 and O.2M from the calculation of 
the observationally accessible quantities for these runs. 

Binary systems were handled by only including the brighter member in the analysis of 
the observationally accessible quantities. This choice is motivated by the fact that for real 
observations, masses will need to be estimated from luminosities. Since binaries in GCs are 
typically not resolved (the separations are below a few astronomical units for the range of 
binding energies considered), we observe mainly single sources. Because the luminosity of a 
main sequence star is highly sensitive to its mass (L ~ M 3 5 ), the lighter star contributes 
very little to the overall luminosity in many cases, and thus the total luminosity will be very 
similar to that of the brighter member. 

To quantify the effects of mass segregation, we examined the radial variations in average 
stellar mass — or equivalently — in the slope of the mass function (if the mass function is 
a power law in the mass range considered, then there is a one-to-one relation between the 
average mass and the slope). As a consequence of energy equipartition, heavier particles 
sink to the center of the cluster within a few relaxation times, increasing the difference in 
average mass between the center and the halo of the cluster. As our main diagnostic of mass 
segregation, we define 



where r hm is the projected half-mass radius of the cluster (computed using only visible stars), 
and (to) is the average mass for main sequence stars with to e [0.2 : 0.8]M Q . Both of these 
measurements are taken from projected radial bins each containing 5% of the cluster's visible 
stars. Because nearly all of the deviation in (m)(r) occurs within this radial range, we are 
maximizing our baseline for measuring mass segregation while using fields with a reasonable 



3.2. Mass segregation 




(6) 
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number of stars. This definition also allows for a straightforward comparison to observational 
data as only two fields per cluster are sufficient. 

Fig. [2] depicts the evolution of A(m) for our iV = 16384 to N = 32768 runs without 
primordial binaries. For each run, we analyze the configuration of the system every 15 
Nbody units (which corresponds to more than 10 measurements per relaxation time). Runs 
with an IMBH are represented as red points, while runs without are blue points. Because 
they were drawn from single-mass King models, our clusters begin out of equilibrium. On 
a relaxation timescale, we see them evolve towards a new quasi-equilibrium state. After 
~ 5t r h(0), most clusters have settled into this equilibrium, with those harboring an IMBH 
showing a smaller amount of mass segregation, i.e. smaller values of A(m). The points from 
the control snapshots 128kk.l, 128kkbh.la and 128kkbh.lb, as well as the two 16k control 
runs (16ks.l and 16ksbh.l) are also plotted in Fig. [2} and are in good agreement with those 
from our models. The data in Fig. [2] come from a variety of initial configurations, not only 
in terms of the particle number but also in terms of initial mass function. The use of a 
differential indicator for mass segregation allows us to cancel out the dependence on the 
global value of (m) (or on the global mass function slope). 

Simulations with a Salpeter IMF (16ks, 16ksbh) contain many more massive remnants 
than the number allowed by a Miller & Scalo IMF (16km, 16kmbh, 32km, 32kmbh). For 
example, a 16k simulation with a Salpeter IMF and m G [0.2 : 1OO]M initially contains ~ 20 
stellar mass black holes, whereas a Miller & Scalo IMF will only have ~ 1. This difference 
causes us to observe a slower growth of A(m) in the Salpeter IMF runs, as a central cluster 
of stellar mass black holes partially quenches mass segregation of visible stars, much like 
an IMBH. However, stell ar mass black holes eject each other from the system within a few 



relaxation times (see also lMerritt et al.l 12004 ). so eventually these systems fully develop the 



amount of mass segregation observed in runs starting from a Miller & Scalo IMF. 

The control runs also reflect this trend. The 128k run without an IMBH (128kk.l) has 
a maximum allowed IMF mass of 15M , meaning there are essentially no massive remnants. 
After 5 relaxation times, it is in good agreement with our simulations drawn from a Miller 
& Scalo IMF, which also produces very few massive remnants. The runs with an IMBH 
(128kkbh.la and 128kkbh.lb), which have maximum initial masses of 3OM , and 1OOM , 
respectively, are both consistent with our other runs with an IMBH. However, 128kkbh.la's 
lower maximum mass results in fewer massive remnants, and thus a value for A(m) closer 
to 16kmbh, which also contains few stellar mass black holes. 128kkbh.lb has more massive 
remnants and behaves similarly to 16ksbh, as we would expect. Finally, we see that our 
16k control runs with a lower IMF cut-off at O.1M fall somewhere between the Miller & 
Scalo runs and the Salpeter runs as far as mass segregation is concerned. Although we draw 
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from a Salpeter IMF in the control runs, the lower minimum IMF mass in these simulations 
results in fewer massive remnants than the other Salpeter runs (~ 8 stellar mass black holes 
instead of ~ 20), but more than a Miller & Scalo IMF would produce. We also note the 
increased scatter in these two runs as a result of excluding the large number of main sequence 
stars with masses O.1M < to < O.2M in the calculation of the observationally accessible 
quantities. 

The situation is very similar when primordial binaries are included (see Fig. [3D : runs with 
and without an IMBH again become well separated after ~ 5t r ^(0). As expected, primordial 
binaries carry lighter particles toward the center of the cluster (e.g., a O.6 + O.2M binary will 
sink to the center like a 0.8M Q single star, but will be observed as a single unresolved source 
with the approximate luminosity of a O.6M star). Hence, mass segregation is partially 
suppressed when compared to the runs where /& = 0. This difference in A (to) is more 
significant in the runs with a Miller & Scalo IMF as compared to runs with a Salpeter IMF. 
Because the runs drawn from a Miller & Scalo IMF lack massive remnants, binary stars 
become more gravitationally dominant, and therefore have a more significant impact on 
the dynamics. Fortunately, the binary-driven quenching of mass segregation is weak when 
compared to IMBH-driven quenching and thus it is possible to discriminate between systems 
with and without an IMBH on the basis of A(m), without the need of assuming a binary 
fraction. 

Combining the data from all our simulations with and without binaries, we can identify 
three regions for the value of A (to) in a collisionally relaxed GC, irrespective of its binary 
fraction: 



• A(to) > O.1M . The system is unlikely to contain a central IMBH. 

• A (to) < O.O7M . The system is a good candidate to harbor an IMBH. 

• O.O7M < A(to) < O.1M . The system may or may not contain an IMBH, depending 
on its binary fraction and on the global IMF (and in particular on the number of 
massive dark remnants). 



In addition, an estimate of the binary fraction based on the presence of a parallel main 



seque nce in the color-magnitude diagram is possible for many observed clusters (IMilone et al. 



20081 ). Application of the mass segregation diagnostic therefore can account for the actual 



number of binaries, resulting in a further reduction in the size of the region of uncertainty. 



Including the set of runs with N = 8192, not shown in the plots but whose A(m) 
is reported in Tab. [TJ we see no trends in A (to) caused by an increase in the number of 
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objects in a cluster up to A^ = 32768. In addition, the iV = 131072 control snapshots are 
consistent with our results, strengthening the independence in the evolution of A(m). An 
increase in the number of particles reduces the deviation from snapshot to snapshot. This 
actually improves the application of this indicator to actual GCs, where the number of stars 
is significantly larger than in our runs. Similarly, we see no significant trends in A(m) caused 
by increasing the IMBH mass up to Mjmbh = 0.03, that is 3% of the entire cluster (see 
Tab. [1]). This suggests that reducing M IMBH below 0.01 would still result in a quenching of 
mass segregation. 



3.3. The origin of IMBH-induced quenching of mass segregation 

The onset of mass segregation along with the initial contraction of the cluster brings the 
most massive stars and remnants into a dense environment. Even in clusters with only single 
stars, the dynamical formation of binaries is inevitable. Because Mjmbh is much larger than 
the typical stellar mass, the IMBH has an extremely high probability of exchanging into a 
binary in a close 3-body encounter. It therefore spends much of its lifetime in a binary or 
stable higher-N system (in more than 90% of our snapshots the IMBH is a member of a 
multiple system). As a result, when massive main sequence stars in our simulations sink 
to the core after energy exchanges with other stars, they are efficiently "heated up" and 
scattered away from the core in encounters with the IMBH and any companions it has. 
The IMBH stochastically moves around the core as a result of these encounters and this 
further enhances the interaction rate because the scatter cone is continuously replenished. 
This mechanism for quenching mass segregation naturally explains the lack of dependence 
of A(m) on the number of particles used and the minimal dependence on Mjmbh, as well as 
suggesting an additional explanation as to why the presence of primordial binaries further 
reduces mass segregation. 



4. Discussion and Conclusions 

We have carried out a large set of direct N-body simulations of star clusters with and 
without an IMBH including a realistic mass spectrum and primordial binaries. While pre- 
vious research has focused its attention mainly on the effects of an IMBH on the surface 
brightness and velocity dispersion profiles of the clusters — signatures that are difficult to 
observe — we searched instead for a different fingerprint of the presence of an IMBH. The 
existence of a massive, central object quenches mass segregation and this effect manifests 
itself in collisionally relaxed clusters through decreased radial variation in the average mass 
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of main sequence stars. This effect does not depend on the mass of the black hole as long as 
it is dominant over the typical mass of a star, nor on the details of the initial configuration 
of the system such as initial mass function, density profile and tidal field strength. The 
amount of mass segregation is only weakly dependent on the binary fraction of the cluster. 
This result allows us to use the amount of mass segregation to separate collisionally relaxed 
clusters with and without an IMBH without the need of additional modeling assumptions. 

A critical requirement for the proposed signature is that the system be well-relaxed, 
so that it has already attained equilibrium with respect to mass segregation. From our 
simulations it turns out that this takes about 5t r h(0). However we can only observe the 
current half-mass relaxation time and this might be shorter than its initial value if the system 
has lost a large fraction of its original mass. To compare our simulations to observations, we 
must thus conservatively restrict ourselves to GCs that: 



1. Are not too influenced by the galactic tidal field (that is, with a tidal to half-light 
radius r t /rhi > 10, which corresponds to tidal fields weaker than the weakest field in 
our simulations). 

2. Have half-mass (3D) relaxation times below ~ 1.5 Gyr, i.e. an age above 8t r h- This 
leaves room for a mass loss of about 50% of the initial mass while still giving an 
integrated age of about 5t r h- In terms of observable quantities this translates into a 
half-light relaxation time below m 1 Gyr. 



Based on the iHarrid (119961 ) catalog, 31 galactic GCs satisfy these stringent requirements 
in terms of relaxation time and r t /r^. The proposed diagnostic could probably be applied 
to more clusters after properly evaluating a dynamical model for their configuration and 
eventually accepting some uncertainty in the selection of likely candidates to harbor an 
IMBH. 

Thanks to the HST treasury survey of galactic GCs, data exist for the cores of many 
clusters that e xplore deep enough to se e main sequence stars down to around 0.2M Q . Along 
the same lines, |Pe Marchi et al.l (120071 ). among others, have also acquired images of clusters 
around the half-light ra dius, in order to calcul ate the global mass function of the system. 
The existing data from |Pe Marchi et al.l (120071 ) are sufficient to apply this diagnostic to a 
few actual clusters and the results from such a comparison will be presented in a companion 
paper. 

In closing, we stress again that while the amount of mass segregation has been proven 
here to be a viable indicator for the presence of an IMBH in simulated star clusters, we 
cannot use this method alone to claim the detection of an IMBH. However, by combining 
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the measure of mass segregation with all other constraints from the velocity dispersion and 
surface brightness profiles, we can select the clusters that seem most likely to harbor an 
IMBH, while at the same time excluding some others from further scrutiny. Once we have 
identified those clusters that are most promising, future observations, such as proper motion 
studies, can focus their efforts to secure a robust detection. 

We thank Enrico Vesperini for useful discussions and suggestions, Holger Baumgardt 
for sharing some of his data with us and the referee for a careful reading of the manuscript 
and for constructive suggestions. This work was partially supported by NASA grant HST- 
AR11284. 
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Table 1: Summary of the N-body simulations. 



Name 


N 


Wo 


IMF 


MjmBh/ Mfot 


M IM bh/M q 


fb 


(A(m)) 


A ( m ) mln 


A ( m )max 


8ks 


8192 


7.0 


Sal 


N/A 


N/A 





0.07 


0.040 


0.112 


8km 


8192 


7.0 


M&S 


N/A 


AT / A 

N/A 





0.13 


0.095 


0.167 


8kbs5 


8192 


5.0 


Sal 


N/A 


N/A 


0.1 


0.09 


0.029 


0.138 


8kbs 


8192 


7.0 


Sal 


N/A 


N/A 


0.1 


0.10 


0.071 


0.143 


8kbsll 


8192 


11.0 


Sal 


N/A 


N/A 


0.1 


0.09 


0.037 


0.116 


8kbm 


8192 


7.0 


M&S 


N/A 


N/A 


0.1 


0.09 


0.048 


0.130 


16ks 


16384 


7.0 


Sal 


N/A 


N/A 





0.11 


0.071 


0.158 


16ks.l 


16384 


7.0 


Sal 


N/A 


N/A 





0.14 


0.112 


0.191 


16km 


16384 


7.0 


M&S 


N/A 


N/A 





0.14 


0.112 


0.174 


16kbs 


16384 


7.0 


Sal 


N/A 


N/A 


0.1 


0.09 


0.060 


0.140 


16kbm 


16384 


7.0 


M&S 


N/A 


N/A 


0.1 


0.10 


0.067 


0.127 


32km 


32768 


7.0 


M&S 


N/A 


N/A 





0.14 


0.108 


0.161 


128kk.l 


131072 


7.0 


Kroupa 


N/A 


N/A 





0.13* 


N/A 


N/A 


8ksBH 


8193 


7.0 


Sal 


0.03 


104.0 





0.05 


0.014 


0.091 


okmbn 


o 1 no 
8193 


/ .0 


M&b 


Am 
0.U1 


jo i 


A 




A A7 
V.VI 


A no/; 

0.036 


A 111 
0.121 


8kmBH 


8193 


7.0 


M&S 


0.03 


129.3 





0.06 


0.011 


0.097 


8kbsBH 


8193 


7.0 


Sal 


0.03 


114.4 


0.1 


0.04 


-0.014 


0.080 


8kbmbh 


8193 


7.0 


M&S 


0.015 


69.5 


0.1 


0.04 


-0.021 


0.079 


16ksbh 


16385 


7.0 


Sal 


0.015 


103.1 





0.05 


0.023 


0.090 


16ksbh.l 


16385 


7.0 


Sal 


0.015 


60.9 





0.06 


0.013 


0.118 


16kmbh 


16385 


7.0 


M&S 


0.015 


128.2 





0.08 


0.027 


0.113 


16kbsbh 


16385 


7.0 


Sal 


0.01 


113.4 


0.1 


0.04 


0.015 


0.078 


16kbmbh 


16385 


7.0 


M&S 


0.01 


141.0 


0.1 


0.05 


0.015 


0.084 


32kmbh 


32769 


7.0 


M&S 


0.01 


240.0 





0.07 


0.051 


0.101 


128kkbh.la 


131072 


7.0 


Kroupa 


0.013 


1000.0 





0.09* 


N/A 


N/A 


128kkbh.lb 


131072 


7.0 


Kroupa 


0.009 


1000.0 





0.06* 


N/A 


N/A 



Note. — We calculated the average, maximum and minimum values for A(m) (in solar mass units) between 
5 and 12 relaxation times for each run. Starred values are not averages, but are from a single snapshot. The 
name of each run indicates: 

• The number of stars in the simulation (8k, 16k, 32k or 128k) 

• The presence of primordial binaries (b if fb > 0) 

• The IMF (m for Miller & Scalo, s for Salpeter, and k for Kroupa) 

• The presence of an IMBH (absent for none, bh for small BH mass, and BH for larger - see also fifth 
column) 

• The value of Wo if different from 7.0 

• Control run with IMF lower cut-off at 0.1M Q (.1 suffix). 
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Fig. 1. — Evolution of the three dimensional half mass radius (r^) and of the core to half 
mass radius ratio (r c /rh) in NBODY units for our N = 32769 simulations with (red curves) 
and without an IMBH (blue curves). The presence of an IMBH prevents core collapse. We 
have smoothed the curves by applying a triangular smoothing window of size 1.0t r / l (0). 
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Fig. 2. — Evolution of mass segregation (via A(m), expressed in M ) across the span of all 
N-body simulations with N > 16384 and /& = 0. Red points are from simulations with an 
IMBH, while blue points represent runs with no massive central object. The runs have no 
primordial mass segregation (A(m) = 0), but on a relaxation timescale, the systems settle to 
a quasi-equilibrium configuration with varying degrees of mass segregation. A central IMBH 
quenches the mass segregation and keeps A(m) < 0.09M o . 
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Fig. 3. — Evolution of mass segregation as in Fig. 1, but for our N > 16384 simulations with 
primordial binaries. Qualitatively, we see that the results are similar to those of the runs 
with single stars, but the equilibrium values of A(m) are marginally lower at later times 
when compared to those where /& = 0. 



